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ABSTRACT 


This  technical  report  summarizes  the  results  of  the  research  we 
have  been  engaged  In  regarding  the  use  of  surface  acoustic  wave  devices  in 
direct  sequence  spread  spectrum  receivers.  The  heart  of  this  research  has 
been  the  use  of  the  device  as  a  real-time  Fourier  transformer.  A  system 
of  this  type  is  sometimes  referred  to  as  a  compressive  receiver,  and  our 
use  of  the  system  has  been  primarily  as  a  means  to  implement  a  narrowband 
interference  rejection  filter.  In  addition,  we  have  studied  many  other 
topics  such  as  rapid  acquisition,  Hilbert  transform  generation,  etc.,  and 
these  topics  are  all  overviewed  in  this  report. 


1.  INTRODUCTION 


During  Che  course  of  this  research,  the  main  emphasis  has  been  on 
narrowband  Interference  rejection  In  direct  sequence  (SD)  spread  spectrum 
communications.  The  research  was  both  analytical  and  experimental  In  nature, 
and  comparisons  between  the  two  sets  of  results  always^ Indicated  excellent 
agreement.  In  addition,  various  topics  other  than  narrowband  interference 
rejection  were  studied,  and  these  will  be  summarized  below.  The  details  of 
these  research  findings  have  been  reported  in  the  form  of  papers  published 
and  papers  presented  in  various  national  and  international  meetings.  These 
are  listed  in  the  appendix.  In  this  report  the  results  are  discussed  briefly 
The  heart  of  virtually  all  the  research  was  the  use  of  a  surface 
acoustic  wave  (SAW)  device  to  perform  a  real-time  Fourier  transform.  The 
output  of  this  Fourier  transformer  was  then  filtered  by  multiplying  its 
output, in  real-time,  with  a  time  waveform  which  represented  the  transfer 
function  of  an  appropriate  filter  (e.g.,  a  notch  filter),  and  this  product 
was  then  inverse  transformed.  Finally,  the  output  of  the  inverse  Fourier 
transformed  was  passed  through  a  conventional  DS  receiver.  This  entire 
operation  was  referred  to  as  transform  domain  processing  (TDP) 

For  the  case  of  narrowband  interference  rejection,  the  interference 
was  taken  to  be  a  single  tone  at  an  arbitrary  frequency  with  respect  to 
the  carrier  frequency  of  the  desired  signal.  Both  experimental  and 
analytical  results  were  obtained  for  the  average  probability  of  error  of 
the  system  when  the  inter ferer  was  fixed  at  a  known  location.  Results  of 
this  type  were  obtained  for  both  rectangular  and  cosine-shaped  window 
functions  at  the  input  to  the  Fourier  transformer. 

In  addition,  it  was  experimentally  verified  that  when  the  location 
of  the  interferer  was  unknown  (indeed,  when  the  frequency  of  the  lnterferer 
was  swept  slowly  across  the  band) ,  an  adaptive  version  of  the  system 
described  above  could  be  implemented  to  provide  close  to  the  same  degree 


2 

of  interference  rejection  as  when  the  location  of  the  interference  was 
known  exactly. 

Using  the  same  type  of  system,  the  effect  of  wideband  Gaussian  inter¬ 
ference  was  analytically  studied.  For  this  case,  the  notch  filter  used  to 
reject  the  narrowband  interference  was  replaced  with  a  prewhitening  filter 
(i.e.,  a  filter  whose  transfer  functions  was  the  inverse  of  the  power 
spectral  density  of  the  noise) . 

As  indicated  above,  topics  other  than  interference  rejections  were 
also  studied.  The  problem  of  rapid  acquisition  of  the  correct  phase  of  the 
spreading  sequence  of  a  DS  system  was  examined,  and  a  new  acquisition  scheme 
using  parallel  SAW  convolvers  was  proposed  and  analyzed.  It  was  shown 
that  a  significant  decrease  in  acquisition  time  could  be  achieved  by  using 
such  a  technique. 

Finally,  the  use  of  TDP  to  generate  waveforms  other  than  Fourier 
transforms  was  considered,  and  techniques  were  developed  to  generate,  in 
real-time,  Fresnel  transforms,  Hilbert  transforms,  and  single-sideband 
signals. 

2 .  NARROWBAND  INTERFERENCE  REJECTION 

Because  the  key  result  of  our  research  is  the  narrowband  interference 
rejection,  in  this  section  we  will  describe  the  operation  of  the  system 
when  used  as  a  narrowband  interference  rejection  filter  and  derive  its 
performance.  All  the  technical  details  of  this  section  can  be  found  in 
[1]  and  [2]. 

The  receiver  to  be  analyzed  is  shown  in  Fig.  1.  The  input  consists 


of  the  sum  of  the  transmitted  signal  ±  s(t),  the  additive  thermal  noise 
n(t)  and  the  interference  I(t).  The  Fourier  transform  of  the  input  is  taken 


the  transform  is  multiplied  by  the  transfer  function  of  some  appropriate 


filter  Hc(w),  the  inverse  transform  of  the  product  is  taken  and  the  result¬ 
ing  waveform  is  put  through  a  detection  filter  matched  to  s(t). 

\ 

Since  the  system  is  linear,  the  three  components  of  the  input  can  be 
treated  separately.  Denoting  any  of  them  by  f(t),  assumed  nonzero  for 
te[0,T],  the  signal  at  point  in  Fig.  1  is  given  by 
T 

f^(t)  *  |  f (t)  cos(waT  -  St2)  cos(wa(t  -  x)  +  B(t  -  x)2)dx 
0 

Simplifying,  this  yields 


T 

fjU)  -  (1/2)  cos(wat  +  St2)  |  f (t)  cos  2Stxdt 

0 

T 

+  (1/2)  sin(wat  +  St2)  |  f(T)  sin  2StTdt 

0 


T 

+  1/2  |  f (t)  cos(2wax  -  2Sx2  +  2Stx  -  wflt  -  St2)dx 
0 

-  (1/2)  F_(2St)  cos(w  t  +  St2) 

R  a 

-  (1/2)  Fj(2St)  sin(wat  +  St2) 


(1) 


(2) 


where  FR(w)  and  F^(w)  are  the  real  and  imaginary  parts,  respectively,  of 
the  transform  of  f(t),  and  the  approximation  used  in  going  from  (1)  to  (2) 
was  to  ignore  the  third  term  of  (1)  which  is  a  double  frequency  term. 
Alternately,  if  it  is  desired  to  exactly  cancel  the  third  term,  one  can 


implement  the  system  shown  in  Fig.  2.  Note  that  this  requires  twice 
as  much  equipment  and*  as  a  practical  matter,  it  is  usually  not  needed. 

Having  transformed  f(t),  it  is  now  desired  to  filter  it  with  a  filter 
whose  transfer  function  is  H^(w).  In  Fig.  1  it  can  be  seen  that  the  output 
of  the  first  chirp  filter  is  multiplied  by  Hc(2Bt).  This  function  is 
related  to  the  desired  transfer  function  H^(w)  *  HR(w)  +  3  H^(w)  hy 

Hc(2Bt)  -  4HR(2Bt)  cos  2Bt  ^  +  4^(200  sin  20t  ^  (3) 

The  terms  cos  2BtT1  and  sin  2BtT1  shift  the  region  the  input  signal  is 
located  in  from  [0,T]  to  [T^^  +  T] .  This  is  necessary  because  the  inverse 
transform  filter  can  be  shown  to  yield  an  accurate  inverse  transform  only 
in  the  range  tetT.,^  +  T]  (see  [1]). 

For  our  purposes,  the  only  filters  H^(w)  that  will  be  considered 
will  be  purely  real,  so  that  Hc(2Bt)  will  be  taken  to  be 

Hc(2Bt)  -  4HR(2Bt)  cos  2BtTx 
resulting  in  a  waveform  at  point  @  of  Fig.  1  given  by 

f2(2Bt)  -  HR(2Bt)  (FR(23t)  cos l (wa  -  28^) t  +  Bt2] 

-  FI(2Bt)  sin[(wfl  -  2BT1)t  +  Bt2]}  (4) 

assuming  the  sum  frequency  terms  at  w&  +  2BT^  are  filtered  out. 

It  was  shown  in  [1]  that  the  region  in  frequency  where  the  transform 
is  accurate  is  the  range  coe[28T,28T^] .  Therefore,  limiting  the  input  to 
the  inverse  transform  chirp  filter  to  te[T,T^]  yields  a  waveform  at  point  <D 


of  Fig.  1  described  by 
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f3(20t) 


T, 

I1 

T  ( 


H  (28x)F  (20x)  cos(wrT  +  &X  ) 


-  Hr(20t)Fi(28t)  sin(wrT  +  B*T) 


.  COs[w  (t  -  T)  -  (t^  -  2tT  +  T^)]dT 

O 


where  w  -  w  -  20  T.  . 
r  a  1 


Expanding  the  terms  in  (5)  and  noting  that  the  only  term  of  interest  for  a 
real  function  f(t)  is  the  term  modulating  cos(wat  —  Bt  )»  the  output  of  the 
lowpass  filter  at  point  of  Fig.  1  is  given  by 

fTi 

f4(t)  =  Hr(20t)  [Fr(2$t)  cos(20  (t  -  T^x) 

T 

-  F^Bx)  sin(2B(t  -  T1)T)]dr  (6) 

Finally,  the  output  of  the  matched  filter  at  point  (i)  is  given  by 

Ti+T 

£<,0^  +  T)  -  f4(t)s(t  -  Tx)dt 


T  T, 


0  T 


|  s ( t)  j  Hr(2Bt)  f(X)  cos  20xXdX  cos  2BtT 


(7  1  \ 

+  Hr(20t)  f (X)  sin  20xXdX  sin  2BtT  V  dxdt 

Lo  J  * 


Letting  w  -  20T  and  interchanging  the  order  of  integrations  yields  (to 
within  the  constant  factor  20) 


T  20T,  T 


+  T)  -  j  |  |  s(t)  [E^CwJfCX)  cos  wX 


cos  wt 


0  2BT  0 


+  HR(w)f(X)  sin  wX  sin  wt]dXdwdt 
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where  SR(w)  +  j  S^(w)  is  the  Fourier  transform  of  s(t).  Assuming  that 
HR(w)  is  bandlimited  to  we[2BT,2BT^]  or  less  yields  (to  within  a  constant 
factor) 

T 

f^  +  T)  -  fCX)^1  {HR(w)S(w)}dX 
0 

T 

» 

-  f (X)  fhR(X)  *  s(X)]dX  (9a) 

0 

where  h^t)  £  JF"1  (H^w))  . 

is  zero  outside  of  we[2BT,2BT^]  is  not 
restrictive,  given  the  fact  that  the  original  transform  was  accurate  only 
in  that  range.  However,  the  significance  of  the  assumption  is  as  follows: 
The  system  output  as  given  by  (9a)  can  be  shown  to  equal 

T 

|  s(t)[f(t)  *hR(t)]dt  (9b) 

0 

This  can  be  seen  by  noting  that  (9b)  is  just  the  output  of  a  filter  matched 
to  f(t)  *hR(t)  when  s(t)  is  the  input  and  when  the  filter  is  sampled  at 


The  assumption  that  H^w) 


t  -  T.  Consequently, 


|  s(t)[f(t)  *hR(t)]dt  -  ~  j  S(w)F*(w)HR(w)e  ^wTe^Wtdw 

0  —00 

S  (w)  F*  (w)  HR  (w)  dw 

—oo 

Using  Che  same  reasoning,  (9a)  can  be  shown  Co  equal 
00 

F(w)S*(w)HR(w)dw 

mmOO 

and  so  (9a)  and  (9b)  are  complex  conjugaCes.  (Recall  from  (3)  chac  H^v 
is  real.)  Since  chey  are  boch  real,  Chey  are  clearly  equal.  Since  (9b;  .j 
Che  ouCpuC  of  a  convenCional  sysCem  (i.e.,  one  in  which  Che  filCering  i. 
performed  direcCly  in  Che  "Cime  domain"),  ic  is  seen  ChaC  Che  syscem  of 
Figure  1  will  perform  very  closely  Co  whacever  parCicular  convenCional 
syscem  ic  is  Crying  Co  emulaCe.  Noce  ChaC  even  if  Che  Cransforms  were 
exacC,  Che  cwo  sysCems  would  noc  perform  idencically  Che  same,  because  Che 
above  analysis  assumes  a  time -limited  input,  and  while  this  is  a  reasonable 
assumption  for  Che  signal  component  of  Che  received  waveform,  Che  noise  in 
most  convenCional  systems  is  assumed  Co  exist  for  all  cime.  Also,  note 
Chat  intersymbol  interference  has  been  neglected.  This  was  done  neither 
as  an  approximation  nor  as  an  oversight.  As  described  in  [1],  to  process 
a  contiguous-time  data  scream,  finice  segments  of  the  incoming  waveform 
have  to  be  cime  gated  and  processed  separately.  If  this  Cime  gating  is 
performed  one  symbol  at  a  cime,  adjacent  symbols  can  be  received  without 
interfering  with  one  another.  The  key  difference  between  the  two  systems, 
however,  is  the  fact  that  components  such  as  ideal  rectangular  filters 
cannot  be  implemented  in  a  conventional  sense,  but  using  the  transform 


t-T 


technique,  filtering  with  an  ideal  bandlimi ted  filter  corresponds  simply 
to  multiplying  the  transform  output  with  an  appropriate  timing  gate. 
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Finally,  it  can  now  clearly  be  seen  how  large  the  range  we [2BT,2BT^]  has  to 
be  relative  to  a  given  signal  and  interference  combination,  since  whatever 
bandwidth  one  would  want  for  H  (w)  in  a  conventional  system  is  the  range 
of  frequencies  [28T,2BT^]  for  which  one  wants  the  Fourier  transform  of 
Fig.  1  to  be  valid. 

To  use  (9b)  to  determine  the  performance  of  a  wideband  system  being 
interfered  with  by  a  narrowband  signal  (e.g.,  a  jammer),  assume  the  trans¬ 
mitted  signal  is  a  binary  waveform  where  each  information-bearing  symbol 
has  superimposed  upon  it  the  seven-bit,  pseudo-noise  (PN)  sequence  -  1  -  1 
-  1  11  -  1  1,  with  Tc  the  duration  of  each  PN  code  chip.  The  noise  n(t) 
is  additive  white  Gaussian  noise  (AWGN)  with  two-sided  spectral  density 
(h0/2)  and  the  interference  I(t)  is  the  tone  a  cos  [(wQ  +  6wQ)  t  +  0]  where 
a  is  constant,  wQ  is  the  carrier  frequency  of  the  transmitted  signal,  and 
0  is  a  random  phase  uniformly  distributed  between  0  and  2ir. 

Each  symbol  is  assumed  to  be  detected  independently  of  other 
symbols.  Hence  the  input  to  the  forward  Fourier  transformer  of  Fig.  1, 
while  certainly  not  time-limited  to  [0,T],  can  be  considered  to  be  observed 
in  disjoint  T-second  intervals.  Therefore  the  input  to  the  system  for  any 
given  symbol  is  effectively  described  by  the  sum  of  that  T-second  long 
pulse  plus  a  T-second  segment  of  both  the  additive  noise  and  the  interfer¬ 
ing  signal.  With  this  model,  the  output  statistic  g(T^  +  T)  equals 
(from  (9a)) 

T 

+  T)  -  j  [±  s (t)  +  n(t)  +  I(t)][hR(t)  *  s(t)]dt  (10) 

0 


VWV.'v.V '.VvV 
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Denoting  the  signal,  noise,  and  interference  terms  in  (10)  by  S^,n^,  and 


In  respectively,  the  noise  term  nn  is  clearly  a  Gaussian  random  variable 


with  mean  and  variance  given  by 


E(nQ)  -  0 


(ID 


and 


v  ■  <V2) 

0 


[hR(t)  *  s(t)]  dt 


(12) 


respectively,  and  the  average  probability  of  error  of  the  system,  as Coining 
equally  likely  symbols,  is  given  by 
2ir 

/  +  in  \ 

(13) 


4TT  /  _  ,  _ 

ij 

0  \  n0 


de 


where 


$(x)  ■  --  -  - 
/2tT 


e-y2'2  iy 


and  where  Iq  is  an  explicit  function  of  6. 


To  compute  either  the  noise  variance  or  either  of  the  deterministic 
terms  in  (10),  one  needs  to  specify  a  filter  h^t)  and  compute 


-1 


s (t)  *  hR(t)  *  F  (S(w)Hr(w)} 


(14) 


For  a  bandpass  notch  filter  with  upper  cutoff  Wq  +  Aw  >  Wq  +  5wq  and  a 


notch  of  width  25w  at  Wq  +  6wq,  the  filter  transfer  function  is  given  by 
Hr(w>  “  p6w<w  "  wQ)  -  P6w(w  -  wQ  -  5w0)  +  P6w(w  +  wQ) 


-  +  Wq  + 


(15) 


For  the  seven-bit  PN  sequence  input  signal  described  above,  it  is 
easily  shown  that  its  Fourier  transform  S^(w)  is  given  by 


I 


S±(w) 


-2  sin3wT  +  2  sinSwT  -  2sin6wT  +  sin7wT 
c  c  c  c 


1-2cos3wT  +  2cos5wT  -  2cos6wT  +  cos7wT 
.  .  _ c _ C  _ c _ c 

w 

Expressing  s(t)  as  s^t)  cos  Wgt,  where  '{S^(w)}  it  can  be 


shown  that 


s(t)  *  hR(t)  -  tR^tjAw)  -  (1/2)  R^t)]  cos  wQt  -  (1/2)  R3(t)sin  wQt 


where 


R^tjAw)  -  -  i  {Si  Awt  -  2  Si  Aw(t  -  3TC)  +  2  S±  Aw(t  -  5Tc) 

-  2  Si  Aw(t  -  6TC)  +  Si  Aw(t  -  7T,)}  ,  (18) 

R2(t)  «-(1/tt)  {Si(6wQ+5w)t  -  Si(5wQ-  w)t  -  2[Si(5w0+6w)  (t-3Tc)  -  Si(6w0“5w)  (t-3Tc)  ] 
+  2[Si[(5w0+6w)(t-5Tc)]  -  Si(6w0-6w)(t-5Tc)] 

-  2[Si(5w0+6w)(t-6Tc)  -  Si(6w0-6w)(t-6Tc>] 

+  [Si(6w0-HSw)(t-7Tc)  -  Si(6w0-6w)(t-7Tc)]}  (19) 


^(t)  --(l/Tr){Ci(6w0+6w)t  -  Ci(<5v0-6w) t-2[Ci(5w0+tSw)  (t-3Tc)-Ci(6w0-6w)  (t-3T£)] 
+  2[Ci(5w0+6w)(t-5Tc)  -  Ci(6wQ-6w)  (t-5Tc)] 

-  2[Ci(6w0+6w)(t-6Tc)  -  Ci(6wQ-6w) (t-6Tc)] 

+  [Ci(6w0+6w)  (t-7Tc)  -  C1(5wq-5w)  (t-7Tc)  ]}  (20) 

In  (18),  (19)  and  (20), 


Also,  for  negative  x,  the  Cauchy  principal  value  of  Ci(x)  is  used. 

Finally,  the  output  of  the  final  matched  filter  can  be  written 
(neglecting  double  frequency  terms) 

T 

» 

gttj+T)  -  [(1/2)  si(t)  +  (a/2)cos(6w0t+6)+n(t)cos  wot][R1(t;Aw)-(l/2)R2(t)]dt 

0 

T 

+  ~  [y  sin(6wgt+6)-n(t)sin  w  t]R,(t)dt  (21) 

0 

Most  of  the  terms  in  (21)  due  to  the  signal  and  the  interference,  and  most 
of  the  terms  comprising  the  variance  of  the  noise  term  in  (21) ,  can  be 
evaluated  in  closed-form  using  the  integral  relations  of  ([3]  p.  633]). 
However,  the  resulting  expressions  are  exceedingly  lengthy  and  consequently 
the  probability  of  error  as  given  by  (13)  will  be  presented  numerically 
below.  As  a  perspective  on  these  results,  if  the  received  waveform  des¬ 
cribed  above  is  detected  with  just  a  filter  matched  to  the  transmitted 
signal  s(t)  (i.e.,  if  no  attempt  is  made  to  remove  the  interference),  the 
average  probability  of  error  of  the  system  will  be  given  by 

2ir 

*  -if  * 

eMF  27r  J 

where  E  is  the  energy  per  bit  of  the  transmitted  signal  and 

T 

I(T)  ■  (a/2)  J  si(t)  cos(5wgt  +  0)dt 
0 

■  !&'"  {[-2  sin36wgTc  +  2  sin5$WgTc  -  2  sin65wgTc  +  sin76wgTc]cos  0 

+  [1  -  2  cos36wnT  +  2  cos56wnT  -  2  cos66wnT  +  cos76wnT  ]sin  0}  (23) 


There  is  one  aspect  of  this  model  which  requires  elaboration,  namely, 
the  reason  for  the  low  number  of  chips-per-bit  (i.e.,  the  small  processing 
gain)  of  the  system.  A  short  PN  sequence  was  used  because  the  shorter  the 
sequence,  the  more  pessimistic  the  results  will  be.  This  follows  because 
the  system  shown  in  Fig.  1  only  "looks"  at  the  interferer  for  T  seconds,  so 
that  the  tone  jammer  actually  has  a  spectrum  whose  main  lobe  bandwidth  is 
2/T.  Since  the  way  this  system  operates  is  to,  for  example,  notch  out  the 
main  lobe  of  the  jammer,  and  hence  unavoidably  null  out  any  signal  energy 
in  that  bandwidth  as  well,  the  larger  the  ratio  of  T/Tc,  the  smaller  is  the 
percentage  of  signal  power  which  is  lost.  In  other  words,  using  a  very 
short  PN  sequence  is  essentially  a  worst  case  analysis  when  comparing  the 
present  system  to  a  conventional  spread  spectrum  system. 

The  results  of  evaluating  the  probability  of  error  expressions 

are  shown  in  Fig.  3.  There  are  eight  curves  labelled  a  to  d  and  a'  to  d’, 

and  Table  I  identifies  each  of  those  curves.  The  notation  P  (a  ■  x)  means 

e 

Eq.  (13)  evaluated  for  a  »  x.  For  these  curves,  the  upper  cutoff  was 
always  taken  to  be  14  tt/T  (corresponding  to  the  chip  rate  of  the  PN 
sequence)  and  the  notch  itself  was  always  centered  at  the  interferer 
frequency  and  had  a  width  of  ±  2tt/2T. 

Comparing  the  curves,  one  sees  that  the  notch  filter  system  almost 
always  outperforms  the  conventional  matched  filter.  The  one  exception  is 
when  the  interferer's  effect  is  initially  not  very  significant  (e.g.,  a 
near  zero).  This  is  because  the  signal  degradation,  due  to  the  bandlimiting 
is  more  harmful  than  the  degradation  due  to  the  interferer.  Clearly,  the 
performance  of  the  system  using  the  notch  filter  could  be  further  improved 

9 

simply  by  increasing  its  bandwidth.  The  key  point  is  that  for  a  strong 
enough  tone  Interferer,  the  notch  filter  in  cascade  with  the  matched  filter 
will  always  improve  the  system  performance. 
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TABLE  I 

INTERFERENCE  LEVELS  USED  IN  SYSTEM  EVALUATION 


Curve 

a 

b 

c 

d 


Description 

Pe(a  ■  0)(i.e.,  no  lnterferer) 
Pe(a  -  2) 

Pe(a  -  4) 

Pe(a  -  8) 
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Finally,  it  is  possible  to  implement  an  adaptive  version  of  this 
system  as  described  in  [2]  and  illustrated  in  Fig.  4.  Its  operation  can 
be  seen  as  follows:  The  upper  branch  envelope  detects  the  Fourier  trans¬ 
formed  input  and  the  output  of  the  envelope  detector  is  fed  into  a  switch 
controlled  by  a  threshold  device.  The  lower  branch  passes  the  Fourier 
transformed  input  directly  to  the  first  multiplier.  The  switch  on  the 
upper  branch  is  set  so  that  any  time  the  output  of  the  envelope  detector 
exceeds  a  predetermined  level,  the  output  of  the  switch  is  forced  to  zero. 
In  this  manner,  the  adaptive  notch  switch  is  implemented. 

Notice  that  the  remaining  processing  as  shown  on  Fig.  4  is  somewhat 
different  from  that  shown  on  Fig.  1,  since  with  the  system  as  implemented 
in  Fig.  4,  the  matched  filtering  is  accomplished  by  Fourier  transforming  a 
reference  signal  s(-t),  multiplying  this  transform  with  the  output  of  the 
first  multiplier,  and  then  inverse  transforming  the  output  of  the  second 
multiplier.  If  the  two  systems  are  compared,  it  will  be  seen  that  they 
are  equivalent  in  the  sense  that  they  perform  the  same  when  operated  under 
similar  conditions  (i.e.,  either  both  in  a  non-adaptive  mode  or  both  in  an 
adaptive  mode) . 

• 

3.  CONCLUSION 

This  report  was  meant  to  summarize  the  results  we  obtained  regarding 
the  use  of  SAN  devices  in  DS  spread  spectrum  communication  systems.  Uhlle 
many  different  topics  were  studied  during  the  course  of  the  research,  the 
most  important  result  was  the  use  of  the  SAW  devices  as  a  real-time  Fourier 
transformer  (sometimes  referred  to  as  a  compressive  receiver)  to  implement 
narrowband  interference  rejection  filters,  and  some  of  these  results  were 
summarized  in  the  previous  section. 
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With  respect  to  results  that  we  obtained  on  the  other  topics  listed 
in  the  Introduction,  a  complete  bibliography  of  all  the  publications 
resulting  from  this  research  is  included  in  Appendix  A.  Hence  the  inter¬ 
ested  reader  can  find  the  details  to  any  of  our  research  results  in  one  or 
more  of  the  referenced  papers.  The  following  two  papers  have  not  yet  been 
published  and  thus  for  convenience  the  copies  of  these  papers  are  included 
as  Appendix  C  and  D. 

Laurence  B.  Milstein,  John  Gevargiz  and  Pankaj  K.  Das,  "Rapid 
Acquisition  for  Direct  Sequence  Spread  Spectrum  Communications 
Using  Parallel  SAW  Convolvers". 


John  Gevargiz,  Michael  Rosenmann,  P.  Das  and  L.  B.  Milstein, 
"A  Comparison  of  Weighted  and  Non-Weighted  Transform  Domain 
Processing  Systems  for  Narrowband  Interference  Excision". 
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-sin  (w  t  -  8t2) 
& 


Figure  2.  Exact  Cancellation  of  Double  Frequency  Term 
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ABSTRACT 

The  effect  of  amplitude  weighting  of  the  received  signal  in  a 
spread  spectrum  transform  domain  processing  receiver  is  studied.  Both 
experimental  and  theoretical  results  of  probability  error  are  given  for 
direct  sequence  (DS)  systems  using  a  63  chip  spreading  sequence  with  non¬ 
contiguous  data.  The  DS  signal  is  received  in  the  presence  of  narrow-band 
interference  and  additive  white  Gaussian  noise,  and  good  agreement  between 
theory  and  experiment  is  demonstrated.  Results  are  given  for  both 
weighted  and  unweighted  input  signals  at  various  notch  bandwidths. 
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1 .  INTRODUCTION 


The  technique  of  transform  domain  processing  for  excision  of 
narrow-band  interference  in  a  spread  spectrum  communication  receiver  has 
been  analyzed  in  [l]-[2],  and  its  processing  gain  advantages  were  described 
in  [3].  Such  a  receiver  has  been  implemented  using  SAW  chirp  devices,  and 
the  experimental  results  of  the  probability  of  error  were  presented  in  [4] . 
It  was  also  shown  that  the  experimental  and  theoretical  results  agree  very 
well. 

The  results  presented  in  [2]- [4]  were  for  an  unweighted  observation 
period  of  T  seconds.  Therefore,  the  transform  of  the  narrow-band  inter¬ 
ference  has  its  spectrum  as  shown  in  Fig.  1.  To  reduce  the  sidelobes,  the 
time-weighting  function  shown  in  Fig.  2  is  used  to  observe  the  input  wave¬ 
form.  It  can  be  seen  from  Fig.  2  that  the  peak  power  limit  of  the  weighting 
function  is  the  same  as  that  of  the  rectangular  weighting.  The  transform 
of  the  weighted  interference  is  shown  in  Fig.  3,  where  it  can  be  seen  the 
first  side  lobe  has  been  noticeably  reduced. 

In  Section  2,  the  implementation  of  the  transform  domain  processing 
receiver  for  the  excision  of  the  narrow-band  interference  will  be  reviewed. 
In  Section  3,  both  theoretical  and  experimental  results  of  probability  of 
the  error  are  presented  with  both  rectangular  and  cosine-squared  time¬ 
weighting  at  the  input.  The  desired  signal  is  taken  to 'be  a  63  chip  PN 
sequence,  and  results  are  given  for  different  values  of  the  excision 
filter  notch  width.  Finally  the  conclusions  are  presented  in  Section  4. 

2.  THEORY  OF  THE  RECEIVER  WITH  WEIGHTING 

The  transform  domain  processing  receiver  to  be  discussed  is  shown 
in  Fig.  4.  For  the  case  of  rectangular  weighting  the  receiver  has  been 


2 


Si 

l 


discussed  in  [1]  and  [2] .  The  received  signal  r(t)  is  composed  of  a 
direct  sequence  PSK  waveform  s(t).d(t),  where  s(t)  is  the  spreading 
sequence  and  d(t)  is  the  binary  data,  plus  the  additive  white  Gaussian 
noise  n(t)  and  narrow-band  interference  I(t).  The  output  of  the  receiver 
of  Fig.  4  can  be  shown  to  be 


g(T  +  T 


i>-J  » 


(t)  [s ( t)  +  n(t)  +  I(t)]  [hR(t)*S(t)J  dt 


where  h  (t)  and  s(t)  are  respectively  the  notch  filter  and  the  reference 
spreading  sequence.  From  the  above,  the  signal  component  of  the  final 
test  statistic,  the  variance  of  the  noise  component  and  the  interfered 
component  can  be  shown  to  be  given  by 


w(t)  s(t)  [hR(t)  *  s(t)]  dt 


f  2  2 

w^(t)  [hR(t)  *  s(t)]Z  dt 


I 

* 

Iq  -  w(t)  Cos [ (cJq  +  6o>q) t  +  6]  [hR(t)  *  s(t) ]  dt 


respectively.  For  a  seven  chip  PN  sequence  the  terms  in  the  bracket 
(i.e.,  h_(t)  *s(t)  is  given  in  [1].  For  other  code  lengths,  the  expression 
generalizes  in  an  obvious  manner. 


3.  THEORETICAL  AND  EXPERIMENTAL  RESULTS 


The  results  are  given  for  a  signal  to  interference  ratio  (SIR)  of 
-20  dB,  where  the  interference  is  a  14.992460  MHz  single  tone.  To  achieve 


adequate  phase  averaging,  the  single  tone  interference  is  phase  modulated 
with  a  phase  excursion  of  2tts  radian  at  100  Hz.  The  curves  are  paramet¬ 
erized  by  normalized  offset  frequency  6o>q  from  the  carrier  frequency  and  by 
the  normalized  single-sided  notch-width  6co. 

Figs.  5-7  show  curve  of  average  probability  of  error  versus  energy- 
per-bit-to-noise  spectral  density  ratio.  These  three  figures  are  given 
for  a  constant  two  sided  notch  width  of  3tt  (90  KHz)  centered  at  14.992460 
MHz . 

Fig.  5  shows  both  the  theoretical  and  experimental  results  of  prob¬ 
ability  of  error  when  the  input  signal  is  time-weighted.  It  can  be  seen 
that  the  experimental  results  are  in  agreement  with  theory  to  within  a 
fraction  of  a  dB.  In  Figures  6  and  7  respectively,  the  theoretical  and 
experimental  results  of  Pfi  for  weighted  and  unweighted  input  signal  are 
presented.  It  can  be  seen  again  that  the  experimental  and  theoretical 
results  agree  very  well. 

In  Fig.  8  the  theoretical  results  of  probability  of  error  versus 
the  notch  bandwidth  for  the  weighted  and  unweighted  input  signal  at  a 
constant  signal-to-noise  ratio  of  10  dB  is  shown.  From  Figs.  8  and  9  the 
optimum  two  sided  notch  width  for  the  weighted  signal  is  4tt  (60  KHz) 
whereas  it  is  12it  (360  KHz)  for  the  unweighted  input  signal.  Also  from 
Fig.  9  it  can  be  seen  that  the  experimental  curve  falls  between  the 
10  dB  and  9  dB  theoretical  curves. 

The  spread  spectrum  code  used  for  this  experiment  is  a  63  chip  PH 
sequence  with  a  chip  rate  of  1.875  MHz.  The  chirp  devices  which  are  used 
to  perform  the  Fourier  transforms  in  Fig.  4  have  a  center  frequency  of 
15  MHz,  bandwidth  of  7  MHz,  interaction  time  of  117  ysec  and  chirp  rate 
of  3.0E  10  Hz/sec.  Finally  the  observation  period  T  is  33.6  usee.  The 


transform  variable  (gj)  is  related  to  time  through  co  =  2Bt  where  B  is  the 
chirp  rate.  Thus  in  the  a)  domain  1  ysec  corresponds  to  60  KHz. 

4 .  CONCLUSION 

Both  theoretical  and  experimental  results  of  probability  of  error 
were  presented  for  a  transform  domain  processing  receiver  using  time¬ 
weighting  of  the  input  signal.  It  was  seen  that  the  theoretical  results 
were  in  good  agreement  with  the  experimental  results  for  both  the  weighted 
and  unweighted  input  signal.  The  most  important  conclusion  is  that  the 
use  of  weighting  does  not  appear  to  buy  anything  in  terms  of  improved 
system  performance. 
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Figure  2 

Trace  A)  Rectangular  Weighting 
B)  Weighting  Function 
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Figure  3 

Transform  of  Weighted  Interference 
(same  scale  as  Fig.  1) 


Figure  7  Probability  of  Error  Versus  E/ri 
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Figure  9 
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ABSTRACT 


In  this  paper,  a  technique  will  be  described  which  uses  multiple 
surface  acoustic  wave  (SAW)  devices  in  parallel  to  reduce  the  acquisition 


time  of  a  direct  sequence  spread  spectrum  communication  system.  Analysis 
of  system  performance  in  both  the  search  and  lock  modes  will  be  presented, 
and  key  quantities  such  as  probability  of  false  alarm,  probability  of 
correct  detection,  mean  dwell  time  and  mean  time  to  lose  lock  will  be 
derived. 
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1.  INTRODUCTION 


In  this  paper,  a  rapid  acquisition  technique  based  upon  the  use  of 
parallel  processing  the  received  waveform  of  a  direct  sequence  spread  spectrum 
communication  system  with  multiple  surface  acoustic  wave  (SAW)  convolvers 
will  be  presented.  Because  the  processing  Is  Indeed  done  In  parallel,  a  de¬ 
crease  In  acquisition  time  from  a  comparable  serial  search  scheme  Is  achieved. 

A  derivation  of  the  probabilities  of  false  alarm  and  correct  detec¬ 
tion  for  both  the  search  mode  and  the  lock  mode  are  presented,  as  well  as  a 
Markov  chain  analysis  of  an  overall  search-lock  strategy.  Expressions  are 
derived  for  such  key  system  parameters  as  probability  of  losing  lock  and 
average  time  to  loss  of  lock.  Finally,  seme  numerical  results  are  presented 
showing  the  variation  of  false  alarm  and  detection  probabilities  with  the 
number  of  parallel  convolvers  and  with  the  ratio  of  received  energy-per-bit- 
to— noise  spectral  density. 

The  paper  is  divided  into  six  sections,  with  a  description  of  the 
system  presented  in  the  next  section,  the  derivations  of  probability  of 
false  alarm  and  probability  of  correct  detection  presented  in  Section  3,  a 
Markov  chain  analysis  of  the  search  and  lock  modes  behavior  given  in  Sec¬ 
tion  4,  and  numerical  results  given  in  Section  5.  Conclusions  are  then 
presented  in  the  last  section. 

2.  SYSTEM  DESCRIPTION 

The  acquisition  scheme  can  be  best  described  by  referring  to  Fig.  1. 
There  are  two  modes  of  operation,  a  search  mode  and  a  lock  mode.  Consider 
the  search  mode  first,  during  which  all  switches  are  in  position  1.  (Ignore, 
for  the  moment,  the  recirculating  delay  line  circuit  at  the  input  to  the 
system.)  Notice  that  there  are  N  convolvers,  each  of  duration  2T  seconds, 
where  it  is  assumed  that M chips  of  the  spreading  sequence  span  T  seconds.  The 
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full  period  of  the  spreading  sequence,  L,  Is  divided  into  subsequences 
each  of  length  M,  and  for  simplicity,  it  is  assumed  that  £  is  an  integer. 

Each  of  the  N  convolvers  has  as  a  reference  input  one  of  the  sub¬ 
sequences  of  length  M.  Initially,  let  us  assume  that  the  total  phase 
uncertainty  at  the  start  of  acquisition,  say  K  chips,  is  spanned  by  the 
MN  phase  positions  of  the  N  convolvers.  (Notice  that  even  if  L  is  much 
greater  than  MN,  this  uncertainty  condition  could  still  be  satisfied.)  This 
guarantees  that  the  correct  phase  position  of  the  received  waveform  is 
"seen"  by  one  of  the  N  convolvers  somewhere  in  the  integration  interval 
of  the  convolvers. 

The  convolver  outputs  are  sampled  at  times  t  *  T  +  jTc,  where  T  *  MTc 
and  j  ■  0,  1,  ...,M-1.  The  length  of  each  convolver  is  chosen  to  be 
2T.  With  these  parameter  values,  if  the  N  reference  inputs  enter  the  con¬ 
volvers  at  t  *  0,  the  unknown  phase  position  will  overlap  the  correct 
phase  reference  of  one  of  the  N  convolvers  at  some  time  during  the  interval 

te[MT  , (2M  -  1)T  ]  ■  [T,  2T  -  T  ].  After  sampling,  the  largest  of  the 
c  c  c 

resulting  MN  samples  is  chosen  as  the  correct  phase  of  the  incoming  wave¬ 
form.  Since  one  of  the  MN  samples  is  the  correct  sample,  it  is  therefore 
possible  to  initially  acquire  in  2T  seconds.  The  actual  mechanism  for 
making  the  final  decision  regarding  acquisition  will  be  described  in  the 
next  section. 

Suppose  now  that  the  initial  phase  uncertainty  K  is  greater  than  MN. 
For  this  case,  the  recirculating  delay  line,  shown  in  Fig.  1,  is  used  at  the 
input  to  the  system  with  both  switches  A  and  B  closed.  On  the  input  side, 
after  the  first  T  second  segment  of  the  received  waveform  is  in  the  system, 
switch  A  is  opened  and  the  T-second  segment  is  recirculated  J  times,  where 
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J  equals  [rrr] ,  and  [x]  Is  Che  smallest  integer  greater  than  or  equal  to  x. 

On  the  reference  side,  the  first  N  subcodes  are  put  through  the  respective 
convolvers  during  the  interval  [0,  T] .  From  [T,  2T] ,  no  reference  input  is 
used,  but  from  [2T,  3T] ,  the  next  N  subcodes  are  used  as  reference  inputs. 

In  general,  the  it^1  set  of  N  subcodes,  i*l,  ...,  J,  is  used  as  reference 
inputs  during  the  interval  [(i-2)T,  (i-l)T], 

Considering  next  the  lock  mode,  switch  A  is  kept  closed,  switch  B 
is  kept  open,  and  all  other  switches  that  were  in  position  1  are  now  put 
in  position  2,  so  that  each  convolver  output  passes  through  the  appropriate 
delay  line.  (Note  that  the  recirculating  delay  line  at  the  convolver  input 
is  not  used.)  In  this  mode  of  operation,  the  N  outputs  are  noncoherently 
summed  to  yield  an  effective  increase  in  integration  time  over  that  used 
in  the  search  mode.  To  reset  the  reference  input  during  the  lock  mode 
simply  requires  separating  successive  subcodes  into  any  of  the  N  convolvers 
by  NT  seconds. 

Notice  that  when  the  system  is  locked  to  the  correct  phase  position, 
all  the  inputs  to  the  delay  lines  have  been  despread  and  hence  the  delay 
lines  themselves  do  not  have  to  be  broadband  devices.  Similarly,  if  an 
incorrect  phase  position  is  being  observed,  while  the  inputs  to  the  delay 
lines  will  no  longer  be  despread,  using  narrowband  delay  lines  will  still 
not  result  in  any  performance  degradation  since  under  this  latter  condition, 
a  signal  component  out  of  the  delay  lines  is  not  wanted. 

3.  PERFORMANCE  ANALYSIS 

th 

Let  the  reference  input  to  the  i  convolver  be  denoted  by 
ra(i^Wt>>  where  P^ct)  denotes  a  segment  of  a  longer  pseudonoise 
sequence  (denoted  PN(t))  from  (k-l)Tc  to  (k-l+M)Tc,  M  is  the  number  of  chips 


in  each  subcode  and  Tc  is  Che  duration  of  each  chip.  Therefore, 

is  a  subcode  in  the  interval  [(i-l)MTc>  iMTc] ,  i  -  1,2,  ...» 
For  simplicity  in  notation  in  what  follows,  the  reference  subcode 
PN^_^M[^(t)  will  be  denoted  by  PN^  (t).  (Note  that  we  have  arbitrarily 
defined  the  start  of  subcode  PN^  ' (t)  to  correspond  to  t*0.) 

To  determine  the  probability  that  the  correct  phase  position  has 
been  found,  assume  that  PN^(t)  is  the  subcode  that  is  currently  being 
received,  and  to  be  specific,  assume  the  received  subcode  PN^ (t)  plus  noise 
enters  the  left  side  of  the  convolvers  of  Figure  1  at  the  same  instant  of 
time  that  the  reference  subcodes  PN^  (T-t) ,  i  ■  1,  ...»  N  enter  the  right 
side.  (Note  that  by  making  the  convolver  length  equal  to  2T,  it  is  not 
necessary  to  have  synchronism  between  the  start  of  a  received  subcode  and 
the  start  of  the  reference  codes;  it  is  done  here  only  for  convenience.) 
Also,  for  simplicity,  assume  J  ■  1.  If  we  concentrate  on  the  ith  convolver 
then  from  [1],  it  is  straightforward  to  show  that  at  any  time  te[T,  2T] , 
the  convolver  output  is  given  by 


U±(t)  -  J  [APNj (x ) cos  wqt  +  n(T)]  PN^R)(-2t  +  ^  +  T  +  t) 


.  cos[w  (x  -  2t)  +  dx  (1) 

o 

where  T^  is  the  length  of  the  convolver,  A  is  a  constant  amplitude,  <{>  is  a 

random  phase  uniformly  distributed  in  [0,  2rr]  and  n(t)  is  AWGN  of  two-sided 
n0 

spectral  density  -y  .  The  specific  limits  of  integration  of  (1)  depend 
upon  the  degree  of  overlap  between  the  received  waveform  and  the  reference 
code.  If  we  assume  the  received  waveform  has  completely  occupied  the  con¬ 
volver  when  the  subcodes  first  enter,  then  for  te[T,  2T] ,  the  limits  are 


0  to  T. 


If  double  frequency  terms  are  ignored  in  the  integrals  and  if  is 
set  equal  to  2T,  then  (1)  reduces  to 


ut(t) 


i 

y  |  A-PH^R)  (-  2t  +  3T  +  T)  PSj(T)  dT  cos  (2uQt:  - 


{  ph<r><- 


2t  +  3T  +  t)  cos  (WgX)n(x)  dx  cos  (2wgt 


{  ph'r)C- 


2t  +  3T  +  x)  sin  (WgX)n(x)  dx  sin  (2wgt  -  <J>) 


Assuming  initially  that  i  ■  j  and  that  the  sampling  time  T  corresponds  to 

s 

precisely  sampling  the  correct  convolver  at  the  peak  of  its  response,  then 
VV  "  [1  AT  +  Nc  (Ts)]  cos  ”  *] 

J 

+  Ns^(Tg)sin  [2w0(Tg)  -  4>]  (3) 

where 


N,  (Tg)  -  j  PNjR^ (x)n(x)  cos  wQxdx 
^  0 


Ng  (Tg)  -  |  PNj(R)(x)n<x)  sin  wQXdx  . 


Both  N  (T  )  and  N  (T  )  are  zero-mean  Gaussian  random  variables  with 
ci  8  si  8 

v 

variance  — j—  .  They  are  also  statistically  independent. 

For  i  +  j  the  output  of  the  convolver  at  the  sampling  instant  is 
given  by 


(6) 


VTS>  -  +  cos£2w0(Ts>  -  ♦! 

+  N  (Tg)  sin[2w0(Ts)  -  <fr] 

where  N  (T  )  and  N  (T  )  have  the  same  statistical  properties  as  N  (T  )  and 
s  s^  S  Cj  s 

N  (T  ),  and  where 
8j  8 

T 

Rij  “  {  PIIiR><t>  PNj(t)  dt  (7> 

0 

Finally,  with  an  appropriate  change  of  notation,  (6)  can  be  seen  to 

apply  to  the  output  of  any  of  the  N  convolvers  at  any  of  the  M-l  incorrect 

till 

sampling  times.  Therefore,  if  the  output  of  the  i  square-law  detector 

at  the  kth  sampling  time  is  denoted  by  S^T  +  kTc) ,  and  if  the  correct 

3T 

sampling  time  (which  equals  y  for  the  specific  timing  being  used  here) 
corresponds  to  k  ■  kc,  then  an  error  will  occur  if 
Si(T  +  kTc)  >  S^T  +  kcTc)  for  any  k,  i  *  j , 

or  if 

Sj (T  +  kTc)  >  Sj (T  +  kcTc)  for  k  #  k£. 

These  conditions  correspond  to  an  incorrectly  matched  convolver  at  any  sampling 

time  having  an  output  larger  than  that  of  the  correct  convolver  at  the  correct 

sampling  time,  and  to  the  correct  convolver  having  a  larger  output  at  some 

time  other  than  the  correct  sampling  time,  respectively. 

Denoting  by  m-l,  . ..  MN  -  1,  the  set  of  samples  of  either 

S^(T  +  kTc)  for  any  k  or  Sj(T  +  kTc)  for  k  +  kc>  and  denoting  S^(T  +  kcTc) 

by  Sq,  then  an  error  occurs  if  Sm  >  Sq  for  any  m.  Hence  the  probability  of 

(s ) 

error  in  the  search  mode,  P  ,  can  be  upper  bounded  by 


■  Pt(S1  >  Sg)  or  (S2  >  Sg)  or  •  or  >  Sg) 

MN-1 

<  l  ?{s±  >  S0) 


P(S^  >  Sj},  then  P^  can  be  shown  to  be  given  by  ([2]) 
j  *  ■j  [1  -  Q(^»  S*)  +  Q(»/a»  */’b)3 


,  f  2a2(A2R.  2  +  A2T2)  -  4a2p. .A2R. .T 


4(p^O 


[4a  (1  -  Pip] 


w 

i»b)  -  | 


2 

a  +  x 


IQ(ax)  dx 


a2  A  n()T 


A  2A  _ 

Pij  "  n0Rij 

and  Iq(x)  is  the  modified  Bessel  function  of  the  first  kind  and  zeroth 
order. 

As  a  special  case,  if  all  the  are  approximately  zero,  then  P 
reduces  to 


P  .1.  2no 

ij  2  * 

and  (8)  reduces  to 


Consider  now  the  lock  mode,  in  which  all  the  switches  in  Figure  1  are 
set  at  position  2.  In  this  latter  mode,  the  N  convolver  outputs  are  all 
appropriately  delayed  and  summed,  but  only  at  a  specific  phase  position 
decided  upon  in  the  search  mode. 

To  analyze  the  performance  of  the  system  in  the  lock  mode,  the 
approximation  that  all  the  *  0  was  made  in  order  to  simplify  the 
analysis.  With  this  approximation,  it  can  be  seen  that  when  an  incorrect 
phase  position  is  being  observed,  the  decision  variable  has  a  central 
chi-square  distribution  with  2N  degrees-of-freedom,  and  when  the  correct 
phase  position  is  being  observed,  the  decision  variable  has  a  non-central 
chi-square  distribution  with  2N  degrees-of-freedom  and  noncentrality 
parameter  aj  -  NA2!2.  Therefore,  the  probabilities  of  false  alarm  and 
correct  detection  in  the  lock  mode  are  given  by 


1  N-l 

■■  ■  y  e 

<2TUT)Nr<N) 


2V 


>4 


respectively,  where  y  is  the  threshold,  T(.)  is  the  gamma-function  and 

I„  ,(.)  is  a  modified  Bessel  function. 

N— 1 

4.  SEARCH  MODE  AND  LOCK  MODE  TRANSITIONS 

In  order  to  complete  both  the  system  description  and  the  performance 
analysis,  the  specific  procedure  for  transitioning  from  one  phase  position 
to  the  next  in  the  search  mode  must  be  specified,  as  well  as  the  procedure 
for  transitioning  from  the  search  mode  to  the  lock  mode.  Towards  this  end, 
a  search-lock  strategy  similar  to  that  described  in  [4]  will  be  used.  In 
particular,  once  a  phase  position  has  been  decided  upon  in  the  search  mode, 
two  more  successive  "hits"  at  that  phase  position  will  be  required  in  order 
to  enter  into  the  lock  mode.  That  is,  once  all  phase  positions  within  the 
uncertainty  range  have  been  examined  and  an  initial  decision  made  as  to 
which  is  the  correct  position,  the  procedure  is  repeated  two  more  times. 

If  the  same  relative  phase  position  is  chosen  both  times,  the  system 
advances  to  the  lock  mode.  Similarly,  two  successive  "misses"  in  the  lock 
mode  will  be  necessary  in  order  for  that  phase  position  to  be  rejected 
and  the  search  mode  reinitiated.  However,  in  the  lock  mode,  only  the 
specific  phase  position  previously  decided  upon  in  the  search  mode  is 
examined,  not  all  possible  phase  positions  falling  within  the  uncertainty 
region. 

A  flow  diagram  of  this  procedure  is  shown  in  Fig.  2.  It  is  well 


known  (see,  e.g.,  [4],  [5])  that  this  type  of  process  can  be  modeled  as  a 


finite  Markov  chain  with  absorbing  boundaries.  From  Fig.  2,  it  is  seen 
that  there  are  six  states  to  this  particular  Markov  chain,  labelled  1  to  6. 
States  1  and  6  are  designated  R  for  "rejection",  meaning  when  either  state 
is  reached,  the  phase  position  currently  being  examined  is  rejected  (i.e., 
the  Markov  chain  has  reached  one  of  its  absorbing  states) .  States  SI  and 
S2  correspond  to  the  two  search  mode  states,  whereas  LI  and  L2  correspond 
to  the  two  lock  mode  states. 

(a)  (a) 

In  Fig.  2,  P  denotes  either  P^v  '  or  1  -  P^  (depending  upon 
whether  the  phase  under  consideration  is  incorrect  or  correct,  respectively) 

(£)  (l) 

and  P^  represents  P^a  or  P^  '  under  similar  circumstances  for  the  lock 
mode.  Using  the  absorbing  Markov  chain  model  for  this  system,  the  canon¬ 
ical  form  for  the  transition  matrix,  described  in  [4]  and  [5] ,  is  shown 


where  I  is  an  identity  matrix,  and  defining  the  new  matrix 

B  -  [I  -  QJ  R, 

it  is  shown  in  [4]  that  the  probability  of  losing  lock  is  given  by  entry 
b2g  of  the  above  B  matrix.  Notice  that  the  notation  b2g  refers  to  elements 
of  the  R  matrix.  That  is,  the  rows  of  B  are  numbered  2,  3,  4  and  5,  and 
the  columns  of  B  are  numbered  1  and  6.  For  the  system  being  considered 
here,  the  probability  of  losing  lock  is  given  by 


26  l-P(l-P) 
s  s 

In  addition  to  the  probability  of  losing  lock,  other  quantities  of 
interest  are  the  mean  dwell  time  in  an  incorrect  phase  position  and  the 
expected  time  to  loss  of  lock,  given  that  the  system  is  in  state  L2  at 
the  correct  phase  position.  From  [4]  and  [5] ,  it  is  shown  that  if  x^  is 
the  mean  time  to  arrive  at  any  absorbing  state  given  an  initial  transient 
state  j  (e.g.,  SI),  then  the  vector  of  elements  x^,  denoted  x,  is  given  by 


[?  -  Q]'1  X  , 


where  X  is  a  vector  whose  elements  are  the  times  needed  to  make  each 
transition. 

Since  the  time  needed  to  move  from  one  state  to  the  next  state  is 
2TK  in  the  search  mode  and  (N+1)T  in  the  lock  mode,  it  can  be  shown  in  a 
straightforward  manner  that  the  mean  dwell  time  is  given  by 

2K(1  -  PT)2  (1  +  P_)  +  (N+l)  P_2(2  -  P,) 


(i  -  pTr  i 


and  the  expected  time  to  loss  of  lock  Is  given  by 


2  “  PL 

T,  -  - — r-  (N+1)T  (24) 

(1  -  PL)2 

5.  NUMERICAL  RESULTS  AND  DISCUSSION 

As  a  specific  example.  It  will  be  assumed  that  a  PN  sequence  of 

length  L*1023  is  being  received  simultaneously  in  N«ll  convolvers,  thereby 

resulting  in  M-93  unknown  phase  positions  in  each  convolver.  If  the 

approximation  leading  to  (14) -(16)  is  used,  namely  that  all  the  p^  of  (13) 

(s) 

are  zero,  the  curve  of  P  w  shown  in  Fig.  3  results.  For  the  same  system 

e 

(A) 

in  the  lock  mode.  Fig.  4  shows  results  of  for  various  values  of  y  (or 

(£) 

equivalently  for  various  values  of  P^a  ). 

The  variations  of  P^a  and  Pd  with  N  are  illustrated  in  Figures  5-9. 
(s) 

Fig.  5  shows  Pfa  versus  E/n^  for  three  values  of  N,  N-5,  10  and  15. 

(S.) 

Similar  results  for  P^  are  presented  in  Fig.  6,  and  Figs.  7-9  show  how 
(I)  (Z) 

Pd  varies  with  P£a  1  for  a  given  value  of  N. 

As  a  perspective  on  the  effect  of  approximating  the  p  by  zero 
for  i  +  j,  the  exact  P^  of  (9)  were  computed,  and  some  of  the  results  are 
shown  in  Fig.  10.  For  the  same  PN  sequence  of  length  L*1023  referred  to 
above,  a  specific  subsequence  of  length  P^93  was  fed  into  all  eleven  con¬ 
volvers  and  the  resulting  1023  values  of  p^  were  computed.  That  is,  the 

(R) 

starting  phase  of  the  reference  input  sequence  PN^  (t)  took  values  0,  93  T 
...  930  Tc,  corresponding  to  (i  -  1,  2,  ...  11),  where  i  is  the  index  of 
the  convolver.  Also,  the  phases  of  PNj(t),  the  input  sequence,  took 
values  0,  Tc,  2  Tc  . ..  92  Tc,  corresponding  to  j  ■  1,  2,  3,  ...  93.  Hence 
was  calculated  for  all  these  combinations  of  i  and  j  over  the  obser¬ 
vation  period  T,  and  these  p^  were  then  used  to  compute  the  corresponding 


values  of  P. .. 


There  ere  four  curves  of  along  with  their  corresponding  values 
for  shown  in  Fig.  10.  These  curves  correspond  to  (i,j)  pairs  of  (2,1), 
(3,1),  (4,1)  and  (1,3).  With  the  notation  introduced  above,  (2,1)  corres¬ 
ponds  to  the  of  the  second  convolver  when  the  starting  phase  of  its 
input  signal  is  zero  and  that  of  its  reference  is  93  Tc<  Similarly,  (1,3) 
corresponds  to  the  first  convolver  when  its  input  has  a  starting  phase  of 
2  Tc  and  its  reference  a  starting  phase  of  zero.  The  maximum  variations  in 
P^j  shown  by  these  curves  at  any  given  SNR  come  very  close  to  the  largest 
variation  in  P^  when  all  possible  (i,j)  combinations  are  considered. 

Also  plotted  on  this  figure  is  the  curve  corresponding  to 

1  -  -j  SNR 

■  0  (i.e.,  ^  e  ).  Upon  comparing  these  curves,  it  can  be  seen 

that,  at  least  for  this  one  example,  the  probability  of  error  generated  by 

approximating  the  by  zero  falls  somewhere  between  the  upper  and  lower 

values  of  P^  for  SNR  between  about  8  and  13  dB,  corresponding  to  a  P^ 

-2  -5 

between  about  2  x  10  and  2  x  10  .  In  particular,  the  approximate  value 

is  within  a  factor  of  two  of  either  the  largest  or  smallest  value  for 
SNR  ■  9  dB,  and  within  a  factor  of  five  (approximately)  at  SNR  ■  13  dB. 

If  the  curves  of  probability  of  error  are  examined,  it  will  be  seen 
that  the  price  one  is  paying  for  the  decrease  in  acquisition  time  is  either 
an  increase  in  probability  of  false  alarm  in  the  search  mode  at  a  given 
E/rig  (relative  to  a  classical  serial  search  strategy) ,  or  an  increase  in 
required  E/Uq  to  yield  a  prespecified  probability  of  false  alarm.  This  is 
to  be  expected,  since  in  the  search  mode,  a  decision  is  made  by  comparing 
JMN  voltage  samples  with  one-another.  That  is,  there  are  JMN-1  incorrect 
phase  positions  present  and  the  probability  that  the  samples  corresponding 


to  at  least  one  of  them  will  exceed  the  value  of  the  sample  corresponding 


to  the  correct  phase  position  increases  as  JMN  increases,  for  a  given 

E/V 

The  advantage  of  the  scheme,  of  course,  is  that  in  21  seconds,  MN 
phase  positions  are  examined.  In  a  standard  serial  search  technique, 
examining  MN  phase  positions  requires  MNT  seconds,  and  so  an  improvement 
of  a  factor  of  MN/2  is  achieved.  However,  it  is  still  necessary  to  have 
an  acceptable  probability  of  error  in  the  search  mode,  and  if  the  E/Hq 
required  to  yield  a  given  level  of  performance  is  too  large,  one  can 
always  increase  the  number  of  states  in  the  search  mode.  While  this  will 

(s) 

not  decrease  P£  on  any  given  try,  it  will  decrease  the  probability  of 
entering  lock  at  an  incorrect  phase  position.  Also,  this  further  illus¬ 
trates  a  point  made  in  [4] ,  namely  that  a  large  integration  time  in  the 
lock  mode  is  desirable  (recall  that  in  the  scheme  described  here,  the 
integration  time  in  the  lock  mode  of  any  one  convolver  is  the  same  as  it  is 
in  the  search  mode,  namely  T  seconds,  but  that  noncoherently  summing  N 
output  samples  has  the  effect  of  yielding  a  larger  integration  time) . 

Since  one  has  to  search  all  the  phase  positions  in  the  uncertainty  region 
in  the  search  mode,  but  one  only  has  to  check  a  single  phase  position  in 
the  lock  mode,  once  one  has  entered  the  lock  mode,  it  is  desirable  to  stay 
in  that  mode  (i.e.,  not  lose  lock)  for  as  long  as  possible.  As  long  as 
the  overall  mechanism  of  entering  the  lock  mode  results  in  a  correct 
decision  with  a  reasonably  high  probability,  the  additional  time  spent  in 
false  lock  caused  by  employing  this  strategy  can  be  kept  to  a  tolerable 
level. 

Finally,  it  is  at  times  more  meaningful  to  consider  such  parameters 
as  meaa-dwell  time  at  an  incorrect  phase  position  and  mean-time  to  loss-of- 
lock  at  the  correct  phase  position  than  it  is  to  consider  the  probabilities 


of  detection  and  false- alarm.  (Which  set  of  parameters  is  more  relevant 

to  a  given  system  is  typically  a  function  of  how  long  the  acquisition 

sequence  is  available,  with  the  former  parameters  being  more  significant 

in  continuously  running  sequences  as  in  certain  navigation  systems 

and  the  latter  parameters  being  more  significant  in  more  bursty  systems 

such  as  packet  radios.)  To  illustrate  the  behavior  of  this  particular 

acquisition  scheme  with  respect  to  mean-dwell  time  and  mean-time  to 

T2 

loss-of-lock,  consider  Figs.  11  and  12.  In  Fig.  11,  —  is  plotted  versus 

the  false  alarm  probability  in  the  search  mode,  with  the  false  alarm 

t4 

probability  in  the  lock  mode  as  a  parameter.  In  Fig.  12,  —  is  plotted 
versus  the  probability  of  correct  detection  in  the  lock  mode. 

6.  CONCLUSION 

A  technique  for  the  rapid  acquisition  of  a  DS  spread  spectrum 
signal  has  been  presented.  The  scheme  relies  on  parallel  processing 
different  subcodes  of  a  longer  spreading  sequence  simultaneously  in  SAW 
convolvers,  and  results  in  a  reduction  of  search  time  which  is  proportional 
to  M  times  N,  where  M  is  the  number  of  chips  being  processed  in- each  .... 
convolver  and  N  is  the  number  of  convolvers.  Results  were  presented  for 
such  key  system  parameters  as  probability  of  detection  and  mean-time  to 
loss-of-lock.  While  only  analytical  results  were  considered  in  the  paper, 
it  is  believed  that  the  system  could  be  implemented  with  current  technol¬ 
ogy  in  SAW  devices. 
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Figure  10.  Probability  of  error -for  various  relative  starting  phases  of 
reference  and  input  for  N  -  11  and  M  »  93  in  the  search  mode. 


Mean-time  to  loss-of-lock  versus  the  probability  of  correct 
detection  in  the  lock  mode. 


